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Abstract 

We develop further a recent dynamical replica theory to describe the dynamics of 
the Sherrington-Kirkpatrick spin-glass in terms of closed evolution equations for 
macroscopic order parameters. We show how microscopic memory effects can be 
included in the formalism through the introduction of a dynamic order parameter 
function: the joint spin-field distribution. The resulting formalism describes very 
accurately the relaxation phenomena observed in numerical simulations, including 
the typical overall slowing down of the flow that was missed by the previous simple 
two-parameter theory. The advanced dynamical replica theory is either exact or a 
very good approximation. 
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1 Introduction 



In a previous paper 0] we introduced a theory to describe the Glauber dynamics of 
the Sherrington-Kirkpatrick (SK) 0| spin-glass model in terms of deterministic flow 
equations for two macroscopic state variables: the magnetisation m and the spin- glass 
contribution r to the energy (for a more general discussion of the SK model and of the 
relevant literature we refer to The theory is based on the removal of microscopic 
memory effects: the only 'knowledge' the system is assumed to have of its past is the 
value of the macroscopic state (m, r). In fact any macroscopic dynamical theory for 
the SK model must contain as dynamical variables, either explicitly or implicitly, at 
least the magnetisation m (which is the relevant observable for many typical spin-glass 
remanence phenomena) and the total energy of the system (in order for the theory to 
reproduce the correct equilibrium equations). Since for the SK model the energy per 
spin is a simple function of the macroscopic state vector (m,r), the theory of [U can 
be seen as the simplest two-parameter dynamical theory for the SK model that has the 
properties of being exact for short times (upon choosing appropriate initial conditions) 
and in equilibrium. 

At a technical level, the resulting formalism is a dynamical replica theory, which at 
fixed-points of the macroscopic flow reduces to the standard equilibrium replica theory, 
including replica symmetry breaking (RSB) a la Parisi || if it occurs. This is in contrast 
to an alternative formalism based on path-integral methods (see e.g. 0, H), where it 
is not yet known how to recover the standard equilibrium results in RSB situations. In 
fact, the potential of the present theory to provide the link between equilibrium replica 
theory and the description in terms of correlation- and response functions (once the 
hitherto neglected microscopic memory effects have been incorporated), we regard as 
one of its most interesting features. 

In [H| the actual macroscopic flow equations were calculated only within the replica- 
symmetry (RS) ansatz. The RS version of the theory was quite succesful at predicting 
the flow trajectories in the (m, r) plane, but also exhibited clear deviations in terms 
of the long-term temporal dependence of the macroscopic state variables m and r. 
These are partly due to the elimination from the dynamical equations of microscopic 
memory effects, and partly an artefact of the RS ansatz. The latter cause for deviations 
can be dealt with by allowing for a breaking of replica symmetry, following Parisi's 
RSB scheme Q; although technically non-trivial, this can be seen as a straightforward 
generalisation of the theory (dynamical RSB within the present formalism will be the 
subject of a subsequent paper Q). Here we concentrate on the more subtle question 
of how to incorporate microscopic memory effects, i.e. on how to generalise the ideas 
in Q to a situation where a macroscopic state specifies the details of the underlying 
microscopic states to a much higher degree. We will show that, by considering as the 
appropriate macroscopic dynamical observable the joint spin-field distribution, one can 
indeed follow the steps in O and arrive at a dynamic replica theory which not only 
inherits by construction the exactness of the previous (simple) two-parameter theory 
in the temporal limits t — > and t —* oo, but also describes the simulation experiments 
accurately, as far as our (limited) data allow us to conclude. The philosophy of our 



3 



approach resembles the one proposed by Horner 0, who also derived a closed diffusion 
equation for the evolution of the joint spin-field distribution. However, at the technical 
level of the closure procedure there are important differences. The present theory is 
not only exact for short times (as is but also in equilibrium, in the sense that the 
full RSB equations are recovered. More importantly, it is constructed in such a way 
that it will produce exact dynamic equations if the joint spin- field distribution indeed 
turns out to constitute a closed level of description. 

This paper is organised as follows. First we generalise the theory of |fl] to the case 
of an arbitrary set of macroscopic observables, and derive constraints on the allowed 
choices for such a set, by requiring exactness in specific limits. We then derive a 
diffusion equation for the joint spin-field distribution, which generates a dynamical 
replica theory. The relevant equations are simplified by making the RS ansatz, and 
their predictions are compared to the results of numerical simulations. We close the 
paper with a discussion of our results and their implications. 



2 Dynamics Sz Replicas 

2.1 Definitions and Macroscopic Laws 

The Sherrington-Kirkpatrick (SK) spin-glass model Q consists of N Ising spins crj 6 
{— 1, 1} with infinite-range exchange interactions J if. 

1 1 

Jij = jjJo + -p^ Jzi i (* < 3) CO 



where the quantities Zij, which represent quenched disorder, are drawn independently at 
random from a Gaussian distribution with zero mean and unit variance. The evolution 
in time of the microscopic probability distribution pt(&) is taken to be of the Glauber 
form described by the master equation 



d 



— p t (a) = \pt(FkO-)w k (F k cr) - p t (a)w k (a)] (2) 



k=l 



in which F k is a spin-flip operator F k &(cr) = 3>(oi, . . . , — a k , . . . , cr/v) and the transition 
rates Wk(cr) and the local alignment fields hi(cr) are 

w k(s) = - [1 - CT fc tanh[/3/i fc (<r)]] hi(a) = ^ Jij^j + & (3) 

where (5 is the inverse temperature, which leads asymptotically to the required stan- 
dard Boltzmann equilibrium distribution Poo(°") ~ ex p[~~ /3H(<t)], with the conventional 
Hamiltonian 

H(a) = -J2 OiJijOj -9j2 a i ( 4 ) 

i<j i 



4 



We now turn to the evolution in time of any given set of I macroscopic observ- 
ables fi(cr) = (Qi(<r), . . . , f^(<x)), described by the macroscopic probability distribu- 
tion T't(O) = J2a Pt{ (T )fi ~~ ^(°")]- We insert the master equation (||), and expand 
the result in powers of the 'discrete derivatives' Af(cr) = Q^FiCr)— f2„(o"), which gives 

5*<«)~ £ ^" {pKn)<E-.('W <»»n,} 
+ 1 E ^^{^(n)(E^( CT ) A f( CT ) A n^))n,t} + o(^ 3 A 3 ) (5) 

where we introduced the sub-shell, or conditional, average 



</(*)>n.* 



£gPt(g)^[n-n(ff)] /W 



If the second (diffusion) term, which is 0(N£ 2 A 2 ), vanishes for N — > oo, equation (||) 
acquires the Liouville form, the solution of which describes the deterministic flow 

jn t = MFurynitr)])^ (6) 

i 

Although exact for N — ► oo (provided the diffusion term indeed vanishes), the set (||) 
need not be closed, due to the appearance of pt(cr) in the sub-shell average. 

There are two natural ways for the set @ to close. Firstly, by the argument of 
the subshell average in @ depending on a only through fi(<r) (now pt{&) will simply 
drop out), and secondly by the microscopic dynamics @ allowing for equipartitioning 
solutions (where Pt(&) depends on er only through f2(<r)). In both cases one obtains 
the correct equations for fi t upon simply eliminating Pt(&) from (|6|). 

2.2 Closed Flow Equation for an Order Parameter Function 

Generalising to the present case, we now make the following assumptions: 

1. The observables fi(<x) are self- averaging with respect to the microscopic realisa- 
tion of the disorder at any time. 

2. In evaluating the sub-shell averages we assume equipartitioning of probability 
within the fi-subsheil of the ensemble. 

As a result pt(cr) drops out, and the macroscopic equations @ are replaced by closed 
ones, from which the unpleasant fraction is removed as in Q using the replica identity 
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We now obtain: 

d r<rS[to-to(<r)]?:iM<r)[to(Fi<r)-n(<r)] \ 

dt 1 [ J2 a 5[n-n(<r)] >{z ^ } 

n 

= HE • • • D£ ^ [n^ 1 )-"^ 1 )] II 5 W (8) 

Our second assumption (equipartitioning) is the dangerous one; its impact on the accu- 
rateness of the theory depends critically on the choice made for the observables f2(er). 
If, however, the observables fi(er) indeed obey closed self-averaging dynamic equations, 
our closure procedure will be exact (see our reasoning above). Requiring the theory 
to be exact in two solvable limits (equilibrium and J — > 0, respectively) imposes con- 
straints on the allowed choices for fi(er). Since in equilibrium we have equipartitioning 
of probability in the energy-subshells (with the Hamiltonian @), and since for J — > 
one obtains closed dynamic equations for the magnetisation, our two requirements im- 
ply fi(er) = (m(cr) , H (a) , . . .) (modulo equivalent combinations). For the SK model 
the energy per spin can be written as 

H{a)/N = -I J m 2 (cr) - 9m(cr) - Jr(er) + ^J /N (9) 

with 

1 1 
m (°") = m E ff < r ( <T ) = 777^7 E aiZ v a i ( 10 ) 



so the choice made in leads to the simplest two-parameter theory that meets our 
requirements of exactness in the two solvable limits. Improving upon Q implies in- 
cluding microscopic information beyond (m,r), i.e. adding observables to the set 
fi(er) = (to(<t), r(er)). Addition of any finite number of observables, although tech- 
nically simple, is not expected to give more than just minor corrections. In contrast 
we choose for the set of observables fi(cr) the (infinite dimensional) joint spin- field 
distribution: 

i 

with the local fields (|3j). Our motivation for this choice is the following 

1. The previous two dynamic parameters m(er) and r(er) can be written as integrals 
over D(^,h; cr), so the advanced theory automatically inherits the exactness in 
the two solvable limits t — > oo and J — > 0. 

2. The order parameter function -D(?, h) specifies the underlying states a to a much 
higher degree than (to, r); i.e. microscopic memory is taken into account. 

3. The microscopic equation @ itself is formulated in terms of spins and fields. 



4. The choice ( |il| ) allows for immediate generalisation to models without detailed 
balance and to soft-spin models. 
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To avoid all kinds of technical difficulties we assume that the distribution ( |T[j) is suf- 
ficiently well behaved; we assume that we can evaluate A(?j h) in a number I of field 
arguments hn and take the limit I — ► oo after the limit N — > oo. We then have 2^ 
observables ?(U (<r) = -D(?> <r), with = 1, . . . ,£ and ? = ±1. In order to work out 
equation ((^) we calculate the discrete derivatives A^ M (<r) = -D(?, h^; Ficr)— D(q, h^; a): 



2 7 2 1 3 

+ jV2E^/[ h r^W " jjWiSK-hiW + OiN-*) (12) 

where primes indicate derivatives (in a distributional sense). Since A^(cr) = 0{N~?), 
the diffusion term in (||) could be 0(1). Explicit calculation, however, will show that it 
vanishes as N~ 2 ; at this stage we anticipate that calculation and assume deterministic 
evolution. To suppress notation we write 



£ JdH f(a, H)D(a, H) = (f(a, H)) D 



In order to see clearly for which terms our two closure assumptions will actually be 
operational, we first work out the exact equation (g). We insert ([14) into (||) and 
retain only the leading 0(1) terms: 

^D t (q,h) = i [l+?tanh(/%)] D t (-^,h) - i [l-<rtanh(/%)] AM) 
+^ J DtCf.^^-e-Jo^anh^H))^] - ^^^^(tanh^H)^/^-^^)])^; 



(13) 



with the sub-shell average 









D(s, hp) — J2j K<yj 6 \K ~ h 3 ( <T )] 





The closed dynamic equation (||) is subsequently obtained from ([13]) by elimination of 
Pt(cr) and averaging over the disorder (using identity @): 

^A(?» = ~ [l+?tanh(/3/i)] A (-?,/») " 5 [l-?tanh(/3/i)] AM) 
+7^ I AM) [/i-e-J o (tanh09fO>A] -T^Y.(( z ii tan HPhi(*))5^ j 5[h-h 1 (*)})} D 

on 1 N ^ N i^j 
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+ J ^ -m D<4 I 1 -** tanh(#K(«r))] ^ CT ,^-^(o-)])) 0t }> + 0(aH) (15) 



with 



((f [a; {z kl }])) D = Hm V • • • EC/V; {*«}] II II 5 

(16) 

Finally, it will become clear shortly that the diffusion term in (^) is relatively simple, 
essentially obtained by replacing zf- — ► 1. As a result all complications are concentrated 
in a single term .4, and we find for A" — > oo the relatively simple final expression 

Y t Dtk,h) = i [l+ f tanh(/370] AK ft) - ~ [l-<rtanh(/%)] AM) 



){*«} 



+^ |AM) [/i-e-J (tanh(^fO>Ut] + D t ] + J 2 [l-(atanh(/3iT)) J ^AM) 

(17) 

with 



J 



A] = - Jim — 7 =E((%' tanh (^(°'))^/[ /l -^( ')]))A ( 18 ) 



3 Replica Calculation of the Flow 

We now turn to the evaluation of {{f[cr;{zki]))D- 
3.1 Disorder- and Spin- Averages 

In the familiar fashion for replica calculations we carry out the disorder averages before 
the spin averages. We remove the disorder dependence through the local fields from 
within the constraining delta functions, by inserting 



1 = nn / dH k 5[m-h k (<T a )\ = nn 

a k 



a k 

so that we can write ( |l6| ) as 



dh k dH k Jh?(H?-h k ((T a )) 
2tt 



((f[<T;{z k i}])) D = lim 



dH L dh 



i 



dH n dh\ E---E e ' 

cr' cr" 



y.y h a \H a -e]- l -lp-Y^.Y h a a 



if 



(fW^izu}] e-*^^^^)^ (19) 
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After symmetrisation due to Zij = Zji and after using permutation symmetry with 
respect to site labels, we find that the disorder averages in (^) involve the following 
two integrals, encountered in the flow term (|i~8|), and in the diffusion term, respectively: 



n Dz v & 



-ijLznY (h a a a +h a (T a ) 



Nz- 



12 



'12 



(20) 



) i<j 

1 1 ,2 

(27r)-2e"2 2 , 



with the Gaussian measure Dz = (2ir)~2 e~^ z , and where we only retained the lead- 
ing 0(1) contributions. Applying (|^) to the trivial function f[a; {zki}} = 1 gives a 
normalisation relation, which we can use to avoid having to perform the remaining 
integrals. As a result we immediately find the diffusion term in (^) to be simply 



•^J^™^) [l-(^tanh(/?iT)) Dt ]} 



(21) 



which proves (|l7|), whereas the remaining disorder-induced flow term ( |i~8| ) remains 
non-trivial. A similar calculation shows, upon substitution of ( |T^ ) into (|^), that the 
second term in the macroscopic stochastic equation (|5|) is of order i^N^ 1 . Given our 
assumption that the limit I — > oo can be taken after the limit N —* oo, this proves that 
the evolution of the distribution Dt(s,h) is indeed deterministic on finite time-scales. 

We now introduce the following set of order parameters (with their conjugates) in 
order to a achieve a factorisation over sites of (|l8l), by inserting appropriate integral 
representations of unity (from which all factors 2ir will vanish in the limit n — > 0): 



m, 



A 



<fo/3({°"}) 



N 



N 



N 



The ^-distribution involving Df(q,h) is also written in integral form. Combination of 
the trio (|Tq,19 ,20 ) then leads to a fully site-factorised expression: 

A[s,h;D] = 

iJ 2 lim lim / dmdmdqdqdQdQdRdRdWdWdb g^E^^-f ^E^^^+^/s^ 

N~*oon^0 J 

ffi"*''] ^ • • • ^ tanh(^ 1 1 )i[/ l -ff 2 1 ] ( S v i ^ 



dH^h 1 
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= iJ 2 lim lim / dmdmdqdqdQdQdRdkdWdWdb e N ^+°W 
{{t&n\i{l3Hi)a a ) M {5[h- H^S^^ha) M + (tanh^iyi)^)^ (<5[^--H"i]^ )(Ti <t q )m} 

a 

with 

oujfj, a a p 

~\ j2 S [QocpQap + RapRpa] + log ^ dHdh(M[H, h, a]) a 

a/3 

(in which (f(cr))er = ^~ n J2a 1 ' ' " Z)<t„ /(°")) an d with the effective single-site measure 
M (all vectors now carry replica- indices only): 

{ f[H, h, a]) M = S""<»>^ »■'>■" J " "" (22) 
WL J/ JdHdhY,aM[H,h,a] V 7 

M[ff , h, <t] = exp < ^jm • a — icr ■ qcr — ih ■ Qh — ih ■ Rcr — i^2 W a ph a hp(T a (jp 
{ <*I3 



- D a {a a ,h^)5[h^ - H a ]+i^h a [H a -6 - J m a ] \ 

afi a ) 



By changing the order of the limits TV — ► oo and n — > 0, the remaining integral can 
be evaluated by steepest descent. It is dominated by the extremum of Vl/ which for 
n > 1 defines a global maximum (the 0(1) term in the exponent will drop out due to 
normalisation, as can be checked explicitly by using the above calculation to rewrite 
«1»d). 

3.2 Simplification of the Saddle-Point Problem 

We can make several immediate simplifications. Firstly, variation of VP with respect to 
W a p, Qa/3, Qa/3 and R aj 3 gives saddle-point equations with which to remove all conjugate 
order parameter matrices from our problem: 

W = iQ = ^J 2 q iq = ^J 2 Q iR=J 2 R j 

Secondly, the scaling freedom in the definition of the conjugate parameters D(<^, h^) 
can be used to take the limit I — ► oo: 

J2 D a (a, - J2 &h.D a (a, h^)f(h^) - J ' dH D a (a, H)f(H) (£ - oo) 
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The result of these simplications and of taking the N — > oo limit is the following: 
A[s, h; D] = iJ 2 lim 



n-»0 



^{(tanh^ffi 

a 

with the effective measure (^), in which M and the exponent \P to be extremised are 
now given by 

M[H h a] = e -^- ff -5 j2ff -Q ff -§ j2 ^-<Zfr-*2„^ a ( CTa '^^ 

* = i ^2 [ dH D ( a ' H )D a (cr, H) m arh a + ]-J 2 ^ [q a pQap + R a pRpa] 

aa J a a p 

+ log JdHdh{M[H,h,a])cr 

with the notation 6 = (9, ... ,8). Next we perform the integrations over the conjugate 
fields h, which leads to an effective measure M involving spins and fields only: 

if[H > a])M = IdHJ2aM[H,a] (23) 

i (24) 

* = dH D(a, H)D a (a, H) + i^m a m a + - J 2 ^ [q a pQap + RapRp a ] 

aa a a/3 

- i log det q + log J dH(M[H, a]) a (25) 

In (p5| ) we have neglected irrelevant constants. At this stage it will be convenient 
to calculate the remaining saddle-point equations, by variation of (25). The first of 
these equations, obtained by variation with respect to D(a,H), enables us to write all 
averages with a single replica-index, involving fields and spins, selfconsistently in terms 
of the original distribution D(a,H): 

D{a, H) = {S aj<7 J[H-H a ]) M (26) 

m a = m = (o)d (27) 
Qa/3 = (cr a ap)M (28) 

l 7 p) (29) 



.Jo 
J 2 



W { (H)d -e-J m + iJ 2 mJ2R, 
B I 7 



P 
J 2 



R a p = -4 ^(^cndH-y-e- Jom+iJ 2 {R)a) 1 ]ap) M (30) 
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J 2 Q a( 3 = "o logdet q - 2 % 1 ' )m 31 

We can now write the flow term .A fllq ) of our diffusion equation as 
= - lim V^- 1 )^ 

rt^O ' — ' 
a/9 



{ (t a nh(f3H 1 )a a ) M (Slh-Hifc^Hp-e-Jam+iJ 2 ^*)^ 
+ (<5[^-iIi](5 ?i(Tl a a }M (tanh(/?i7 1 )[^ /3 -^-Jom+iJ 2 (i?tcr) /3 ]) M } (32) 



3.3 Equilibrium 

In equilibrium, we know that the microscopic probability distribution is of the Boltz- 
mann form, Poo(cr) ~ e -P H ( a ) . Therefore, the present constraint restricting micro- 
states under consideration to those with the same joint spin-field distribution, must in 
equilibrium reduce to a constraint selecting states with the same energy. We will now 
make the ansatz 

D a (a,H) = Upa[H + 6] (33) 

and show that it indeed corresponds to a stationary state for our diffusion equation (p]), 
in which one recovers the familiar equations from equilibrium statistical mechanics; i.e. 
the full (RSB) order parameter equations [J as well as the equilibrium local field 
distribution |§. 

We first turn to the saddle-point equations. Given the simple expression ( |33|) we 
can perform the field integrations, with the result 

If 1 
* = ~oP n Yl dH D(a,H)a[H+6] + i^m a m a + - J 2 ^ [q a pQ a p 

a J 

a J a aj3 

+ log(e [/36 ^ '■Jom-irh]-a+±J 2 a '-[^q-Q-ipR^a^ ^ 

(again we forget about irrelevant constants). The remaining saddle-point equations 
become ^ 

™a = -^il3Jom a Q a/3 = --(3 2 q a/ 3 R a/ 3 = -if3qp a (35) 



( aae /3[j m+0]-(T+±(i3j) 2 (T-q(T} cr ^ ge /3[J o m+0]-cr+±(j3J) 2 cr-q(Tj 



, e l3[J rn+0\.<T+\{j3jy<T-qcr\ Hap ^[j m+6\-a+\{i3j) 2 a-qa 



(36) 

which are the familiar equationspl as obtained by an equilibrium (thermodynamic) 
analysis. With the relations (33],35|) we can simplify the effective measure M consider- 
ably: 

if[H > a])M ~ Jdz e -^- lz (e^lJ m + e + Jz ])(T (3?) 



12 



(with m = (m, . . . , m)) This simplified measure obeys useful relations like 

(a a f[H; {a^ a }]) M = (tanh(j3H a )f[H; {a^ a }]) M (38) 

((H a -J m-e)f[{H^ a };a]) M = /3J 2 £ q a p(apf[{H^ a }- <t]) m (39) 

& 

In particular we now find m = (tanh(/3iT))£). If we combine the expression (|37| ) with 
(|26|), sum over the remaining spin variable er and perform the integration over z, we 
are led directly to the equilibrium expression for the local field distribution as obtained 
in [|: 

D(h) = lim / ™ e -hJ^k(h-j om -e)^ 2 V 

v ; n->o./ 2vr , e /3[J rrHO]-er+±(j3J)2cr-q<T\ 

Next we show that the choice (|33|) corresponds to a fixed point of the diffusion 
equation (p7|), i.e. that ^D t {<;,h) = for all {$,h). In the right-hand side of (|i"7|) 
the first two terms trivially cancel, which follows from applying to ( |26|) the identities 

= and d|): 

[l+ ? tanh(/3/i)]D(-^,/i) - [l-?tanh(/?/i)]D(?, fc) = [tanh^/t) -cr a ]) M = 

Equivalently: 

D(q,h) = ~[l+stsaih(ph)]D(h) 



We use ( |35| ) and (|38|,39) to rewrite (|32[). In doing so we will also use equilibrium 
relations like 

^E^L = (tanh(pH){H-J m-9)) D 

a 

which can be derived directly from the equilibrium saddle-point equations (see e.g. ||]). 
The result is: 

A[s,h;D] = -{h-9-J m)D(^,h) - /3J 2 [1- (tanh 2 (/3if)) D ] ? Z%./i) 

+ [l-{tuDh 2 (pH)) D ] lim o E(9 _1 )i7^^-^i]^i[^7- J o^-e]}M (40) 

7 



In order to combine the flow term A in (17) with the diffusion term, we apply (|37|) to 



equation (|26|) and calculate the field derivative: 

9 9 Sdz 5[h-J m-6-Jz a ]^- {e-^-^iS^e^om+e+jz^X 

J 2 —D(q, h) = J lim 1 — ^ >- 

dh n-*o jd z e -±z-q- x z / e po-[j m+0+jz]\ 

= (3J\D(c;,h) - limV(g- 1 )i 7 (<5[/ t - J ff 1 ]5 ?)(T1 [F 7 -Jom-0]}M (41) 

n — 4.(1 ' * 



n->0 ' 

7 



Insertion of (fl0|) and fl4l|) into the right-hand side of (|17| ) leads to the desired result: it 
exactly vanishes. This completes the proof that the standard thermodynamic equilib- 
rium state, as calculated within equilibrium statistical mechanics, defines a fixed-point 
of our diffusion equation (|T7|). Note, however, that this leaves open the possibility of 
existence for stationary states other than the thermodynamic one. 
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4 Replica Symmetric Flow 
4.1 Derivation of the RS Equations 

In order to proceed further in evaluating explicitly the saddle points we now make, as 
a first step, the ergodicity or replica-symmetry ansatz (RS). All order parameters with 
a single replica index are assumed not to depend on this index; all order parameter 
matrices are assumed to have entries which depend only on whether or not they are on 
the diagonal. With a modest amount of foresight we put: 

m a =m q af3 = (l-q)5 a/3 + q 

m a = in R a/ 3 = i(l-q)[R 5 a/ 3 + R\ (42) 

D a {a,H)=ilog X ((T,H) Q af3 = Q 5 a(3 +qR 2 -2(l-q)RR -Q 2 

which implies (q~ 1 ) a (3 = (1 — l)~ 1 [^af3 — ^(l - q)~ 1 ] + 0(n). Working out the RS version 
of the extensive exponent (|25| ) gives: 

Iim — = -mfi- ^log(l-g) - -p— - \j\l-q)Q 2 - J 2 (l-q) 2 [R 2 + 2R R] 
n— »0 n I 2(L — q) 2 

-E l dH D & F ) lo § ^ F ) + lim n - lo S / dH(M RS [H, (43) 
•~ J n— »0 n J 

with 

M RS [H,a] = J] | x(CTa) ^ a) /^-i^^^o-)V^(^Jom) CTa | 

OL ' 

e TJT^^ H ^« m )\ 2 + \ j2 ^ (44) 

We can obtain a factorisation of Mrs [if , cr] with respect to the replica labels by intro- 
ducing appropriate Gaussian integrations: 



with cos(2^) = C ' j2\J AB . Application the above identity to ( f44j) leads to an expression 
for (^3|) in which we can take the remaining limit n — > 0. We use the definition of the 
angle <j> to eliminate the order parameter R from our problem and write the averages 
over the two Gaussian variables x and y as ((. . -)) X y The final result involves an effective 
measure M R s[H, a] without replica indices: 



lim n~~^ = -ra/x- j-log(l-q)- Q , -J 2 (l- 
n-+0 n 2 2(1 — q) 



^Q 2 + {l + q)R 2 ) + 2R Q^q'cos(2(f ) ) 



E JdH D(a,H)log X (<T,H) + ^log J dH (M RS [H , a}} (45) 
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with 

r , , , [Mr {H rfif,° m ] i +Ro(H-0-J m)(T+ {H-8-J m)(x cos sin 0)+JQcr(z cos tp-y sine 

iWRS[-n, trj = x(<7 5 -n Je 2J (1-9) -n 1 <?j 

(46) 

We write averages with respect to this final measure Mrs[H, a], which are parametrised 
by the Gaussian variables x and y, as 

{ f[H,a}), = f dH ^M RS [H,a]f[H,a] 



JdH J2a Mrs [H, a] 

To further reduce our future bookkeeping we derive two useful relations by partial 
integration over the Gaussian variables: 

x(f[H,*])*) = ^^^(f[H 1< r]B) 1t -(f[H, t r]) it (B) 
+ JQcoB<f,((f[H,a]a) it -(f[H,a]) it (a) 

'ixy 

y(f[H,a]),))^=^^^(f[H,a]H),-(f[H,a]UH 
JQ sin 4> ( (f[H, o)o)± - (f[H, 

'Ixy 

Functional differentiation of (|45|) with respect to the function x gives the RS saddle- 
point equation 

D(,,h) = ^(5[h-H]5^)^ (47) 
which implies, as expected, the general relation 



{f[H,a])*) =(f[H,a]) D 

II xy 

Differentiation of ( |45| ) with respect to the parameters {q,m, /i, Rq,Q,4>} and repeated 
usage of the above bookkeeping identities gives the remaining RS saddle-point equa- 
tions: 

m = (a) D (48) 



< = {™L (49) 

2J 2 R (l-q) 2 = ( a (H-J o m-0)) D - ( ( a ),{H- J m-9)A (50) 



w // xy 

2J 2 {l-q) [R (l + q)+Q^cos(2<f>)} = (a(H - J m-8)) D (51) 

H + J R m = — -(H-J m-9) D (52) 
J 2 {l-q) 



15 



J^a-JHI-Q? 



Q 2 +4qR 2 o +8QRoVqcos(20) + ((H -J m-6) 2 ) D 



1+g 
1-4 



({H-J Q m-ef)*-(H-J m- 



(53) 



We now use the RS ansatz to perform the n — > limit in the flow term A (32) of 
our diffusion equation (|T7|). Note that, due to n — > 0, we may deal with averages over 
the original measure M which involve two replica indices (such as those encountered in 
(|32|) in the following way: 



{f [H Q , Va}g[H , a p ]) M - 5 a p(f[H, a]g[H, a]) D + (l-5 af3 ) (( (f[H, a))*(g[H, a]) 



With this identity we can work out (|32"1), We use the short-hand Qy / gcos(2</>) = (1- 
q)R\i and find after some bookkeeping and some re-arranging of terms: 



(l-^) 2 A RS [^h;D] = {2q-l h) 



(h- J m-6){t&nh(PH)a) D +s(t&nh(pH)(H - J m-6)) D 



-qD(s,h) 



+2sJ 2 (l-q) 2 D(<;,h) 



(h-J m-8) (<tanh(/3iT))*(ff)*)) +q ( (tanh(/3if))*(ff- J m 



■ni 



xy 



(R 1 +R ){tanh((3H)a) D -R 1 ( (tanh(pH))*(a 



xy 



+ ({6[h-H}S^U{H-J m- 



xy 



(tanh(£ff))*<(j)*)) -q(tanh(PH)a) D 

xy 



+ ((5[h-H]5^Ua) 



(t&nh{(3H))±(H-J m-8)±)) - q (t&nh{(3 H) (H - J m- 6)) D 

xy 



+2J 2 (1-^) 2 ^{S[h-H]6^{a)A (Rt-Ro) ^(tanh(/?tf )>*(<r)* j> -^(tanh^^o 

(54) 

In RS approximation the evolution of the joint spin-field distribution is described by 
equation (|i~7|), in which the disorder- induced term A is given by (|54|). Evaluation of A, 
in turn, requires solving the set of saddle-point equations (47[53), at each instance of 
time. 



4.2 The AT Instability 

The de Almeida-Thouless (AT) instability || marks the instability of the RS solution of 
the saddle-point equations to the so-called replicon mode. This leads to a second order 
transition away from the RS state to states with broken replica symmetry (RSB). 
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Unlike the standard equilibrium calculations, we here have to worry about replicon 
fluctuations with respect to three replica matrices: 



Qa/3 -> Qa/3 + &<la(3 



$q 



£ ^ = Y • 



laf3 



Qa/3 — ¥ Qa/3 + $Qaf3 



SQaa =0 J2 S Q<*0 = £ 6 Q"0 = ( 55 ) 



Ra/3 RaP + iSRa/3 



SRo 



£ 5R a p = £ SRap = 



with 5q a p = (5^/30, 5Q a f3 = ^Q/3q and <5i? aj g = 8Rp a . As usual the replicon fluctu- 
ations satisfy convenient matrix commutation relations, like [g RS ,<5g] = [Q RS ,5Q] = 
[J2rs, SR] = 0. The AT instability corresponds to a zero eigenvalue in the spectrum of 
the Hessian (i.e. the matrix of second derivatives) of at the RS saddle-point. How- 
ever, since the R a p are conjugate order parameters, acquiring an imaginary value, the 
naive picture of this zero eigenvalue signalling the bifurcation of a local maximum, need 
not be true. We can avoid all such subtleties by following the alternative procedure: 
to consider RSB fluctuations only after elimination of the conjugate order parameters 
R a p with equation (|30|). This is equivalent to first working out the variation in ^ (p5|) 
for the case where all fluctuations ( |56| ) are independent, followed by a projection onto 
the subspace defined by 



Expansion of (25) around the RS saddle point, the first non-trivial order of which 
must by definition be quadratic in the replicon fluctuations, gives 

* - ^rs = \{G 2 ) M + \j 2 £ 5q a p5Q aP - J 2 £ 5R lp ~ j2R o Y SR »P^ap 

~ a/3 aft a(3 



1 



+ «£&& 



ail 



1 



2(1- 



■+3J 2 R 2 - 



1 



(H 2 )*-(H)l 



+ 0(6 3 ) (56) 



with 



G = --J 2 a 
2 



5Q+2RoSR-mSq 



J 2 (l-q) 3 

(H-J m-6) ■ 5q(H-J m-0 



2J 2 (i- q y 



a ■ [5R-R 5g] (H-J m-6) 
1-9 



In order to evaluate the term in (|56D that involves G, we note that in RS saddle-points 
and for indices a ^ (5 and 7 7^ A: 



+ 6 aX s^ (( </fc>*(^>*+</>*<ff>*W*<*>*-(/fc>*^>*<^*-<^>*</>*(fc>. 



■i-y 



1 xy 
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+ terms with less than two <5's 



Only the terms with two Kronecker S's can contribute to (G 2 }m, due to the specific 
properties of the replicon fluctuations. We now obtain 



i(G 2 >M = ^ 4 (([l-Wf 2 



J2[SQa/3+'2Ro5R a f3- R 2 5q a(3 } 2 

X V a/3 



4J 4 (l-g) 4 



X V a/3 



+ 



2(l-g)2 



2(1- 



l-(a)l (H 2 )*-(H)l +[(aHU-(aU(H) ic ] 2 } X^i^-i^g^] 5 



[(aH)+ -(a}*(H}*}} 2 } £ ^ [5Q a/3 +2R 8R aP - R 2 Sq af3 ] 



J 2 



1-q 

+ 



\-{a)l [(aH^-iaUH)*]) Y,i SR ^- R o 6 9 a /3}[5Q a /3+2R dR a(3 -R 2 5q a p] 

ffxy a p 

I'xy Q(3 



J 2 (l-qf 



The various combinations of matrix fluctuations can be somewhat disentangled by 
introducing the transformation 



R 



1 



5Q = -R^Sq-2^Sr+-^5k 6R = R 5q+-5r 

In addition this renders all fluctuations dimensionless. Expression (|56|) now acquires 
the form 

/ Sk a p\ ( 5k a/3 

* " ^RS = J2\ S 9°P I M 
a/3 \5r al3 



Sq a p | + 0(5 3 ) 
in which the entries of the symmetric 3x3 matrix M are 



(57) 



[i-<*>a 



21 2 



M 



12 



Mi 3 = M 31 
M 22 
M 23 = M 32 



1 igvgi] 

1 J2(i=5) 



(58) 
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We now use (f30() to eliminate the conjugate order parameters R a p from our equations. 
In the space of RS saddle-points and replicon fluctuations we satisfy [R, q] = 0, so 
equation ([}(]) simplifies to 

i 

(qR) a/ 3 = -^j 2 -{(H a -J m-9)ap) M 

which after some algebra translates into the following constraint on the replicon fluc- 
tuations 

M 31 5k + M 32 5q + M 33 5r = 

The stability of the RS saddle-point against replicon fluctuations is now controlled by 
a symmetric 2x2 matrix M: 



* - ^RS 




+ 0(5 3 



M 



1 -M31/M33 



M 



1 -M 32 /M- 



33 



1 







1 



-M 31 /M 33 -M 32 /M : 



33 



M n -MyM 33 
M 12 -M 13 M 32 /M 33 



M 12 -M 13 M 32 /M 33 
M 22 -M% 3 /M 33 



(59) 



Due to the curvature sign change of the second derivative of the analytic continuation 
to n — > of the saddle-point that maximises ^ for n > 1, will minimise for n < 1. 
This is emphasised explicitly by the summation over n(n — 1) non-trivial terms (all 
index combinations with a 7^ (3) in (|57|). We can conclude that the AT instability 
occurs when the largest eigenvalue of the matrix M is zero. 



4.3 Equilibrium 

From our previous result, the confirmation that the general (RSB) thermodynamic 
equilibrium state is a stationary state of our flow equation (|T7|), it follows that the 
same must hold within the RS ansatz. We will now show this explicitly, as a non- 
trivial consistency test (rather than a new result). The previous ansatz (33) translates 
into 

x ^ H) = e ^(H+e) (6Q) 

Due to ( |60|) the measure Mrs in (^) becomes a Gaussian function of the fields, which 
enables us to perform the field integrals and work out the RS saddle-point equations 
HBl). The result is: 



M = 2^ J o m 



Q = 



Rn 



(61) 
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m = J Dz tanhf3(J m+9+Jz^/q) (62) 

q = Jdz tanh 2 l3(,hm + e + Jzy/q) (63) 

which are the familiar RS equilibrium saddle-point equations, as first obtained in Q. 

We now turn to the right-hand side of equation Since Q = the original two 

Gaussian variables (x, y) in (|45|) are replaced by a single one, z. With (p0|) and (|35|) we 
can simplify the effective measure Mrs as in the full RSB case, leading to 



(f[a,H]Ug[a,H])A - / Dz (f[a, H\)M°, H}} 



■'II 



_ Y, a eM J ° m+e+Jz yffl fDw f[a,J m + e + Jw^T=q-+Jz^q+pj 2 a(l-q)} 
^ L<7 ' J '*~ 2 cosh /3[J o m+0+J^] 

In particular: 

Psh- i/3 2 J^(l- g )- 57 ^L_ j om -e~ Jz^/q] 2 

^[/i-i?]^ = (64) 

V L J ' ' 2J^2Tr(l-q) cosh (3[J m+9 + J z^q] y ' 

The dependence of ( |64| ) on q only through a factor e^ qh immediately ensures that the 
first two terms in the diffusion equation ( |i~7| ) cancel. Since this happens even before we 
carry out the Gaussian average, we may write 

{5[h-H]5^)+ = i[l+?tanh(/3/0]-D(M 

implying relations like 

(af(H% = <tanh03iO/(iO>* 
The building blocks of ( |5"4| ) thereby become 

(cr)* = tanh/3[Jom + + J#i/(z] 

(H-Jom-9)* = (3J 2 (l-q)(a), + Jz^~q 

J D*(tanh(/3#)>*(*>* = q 
J Dz(tanh(f3H)}*(H-J m-6}+ = 2[3J 2 q{l-q) 

J Dz(tanh((3H)(H-J m-9)}* = f3J 2 (l-q 2 ) 
We will also need the following identity, obtained by partial integration over z: 

J Dz z{5[h-H]5, <a ), = y h)(h-J m-6) - (3J^q(l-q) J Dz{5[h-H]5^),{o)* 

(65) 
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We now have the necessary tools to analyse with minimum effort the complicated terms 
in our diffusion equation, given the ansatz (60). The combined now terms in ([17]) can 
be simplified to 



D{^h){h-J m-e)D^ 1 h)+A KS [^h;D} = 1- (tanh 2 (/?tf)) D 



l-2q 
O^ 2 



D(s,h)(h-J Q m-0) 



-PJ\D^, h) + 



l-Q 



Dz(8[h-H}5^ a )*(a)+ + — ^ J Dz z(8[h-H}S,, a ) 



In order to evaluate the diffusion term in ( |17D we calculate the field derivative of D(s, h), 
using 



J 2 ^-D{q,h) ={3j\D{q,h) 



h-Jom-9 M , -V? 



dh v " ' ' 1 — q v " ' ' 1 _ q 

The full right-hand side of (|T^) can now be written as 

'q(3J 2 



RHS 



l-(tanh i (/3i/)) D 



1- 



Jy/Q 



Dz z{5[h-H]5 q 



r D(q,h)(h-J m-6) 



+ 



J Dz z{5[h-H]5 q 



(1-qf 

(by virtue of the identity (|65|)). The RS equilibrium state obtained in Q therefore 
defines a stationary state of our RS diffusion equation (|17|,^4|) . 

Finally we turn to the AT instability, which we found to occur when the largest 
eigenvalue of the matrix M (^) is zero. We can use the various identities, derived for 
the thermal equilibrium state, to simplify the matrix elements of M considerably: 



M 



cq 



2(/?J) 2 A-l 



( A i[l-(/?J) 2 A] 
i[l-(/?J) 2 A] ±(/3J) 2 [l-(/3J) 2 A] 



\2 



with 



A = Jdz cosh" 4 /3[J m+e+Jz^/q] 



The AT instability, as calculated within equilibrium statistical mechanics |J, occurs at 
(PJ) 2 A = 1. Substitution of this condition into our expression for iW eq immediately 
leads to the desired result: the two eigenvalues of M eq are {—A, 0}, so the two conditions 
for the AT instability coincide. 



5 Comparison with Simulations 

In order to verify the predictions of our theory we here compare the results of solving 
numerically the (macroscopic) diffusion equation (|l7|), in which the disorder-generated 
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term A is calculated within the RS ansatz (|54|), with the results of performing numerical 
simulations of the discretised version of the underlying microscopic stochastic dynamics 
(HH). Solving the diffusion equation (|l~7|), requires making a discretisation not only of 
time, but also of the joint spin-field distribution, i.e. replace the two continuous func- 
tions Dt(±l,h) by two histograms. Furthermore, at each time-step we have to solve 



the RS saddle-point equations (47-53), which involve nested Gaussian integrations. It 
will be clear that the solution of equation ([I?]) requires a significant computational 
effort, even within the RS ansatz, which is reflected in the scope of the experiments 
described in this paper. We restrict ourselves to describing the evolution of the system 
in zero external field (9 = 0), following initial states with individual spin states chosen 
independently at random, given a required initial magnetisation. Following the various 
experimental protocols that show spin-glass ageing phenomena, such as relaxation fol- 
lowing cooling in a small field, and relaxation with intermittant temperature increases 
or decreases, we consider to be beyond the scope of this paper. 

5.1 Transients 

First we study the relaxation of the system on short time-scales. We measure as a 
function of time the magnetisation m, the energy per spin E, and the two distributions 
Dt(±l,h). Note that the full local field distribution Dt(h) is just the sum Dt(l,h) + 
Dt(—1, h). The numerical simulations were carried out with systems of N = 8000 spins, 
following randomly drawn initial states. The results of confronting our theory with 
typical simulation experiments, for relaxations at T = 0, are shown in figures [l] and 
for Jo = (left pictures) and Jq = 1 (right pictures). In figure [j] the top graphs 
represent the magnetisation m and the bottom graphs represent the energy per spin E; 
for the two initial conditions mo = and mo = 0.3. Figure || shows the corresponding 
distributions D(a, h) for one particular choice of initial state (Df(l, h): upper graph in 
t = window, right graph in t > windows; Z)j(— l,/i): lower graph in t = window, 
left graph in t > windows). For Jo = 1 we were not able to calculate the solution of 
equation (|l~7|) up to t = 6, due to the critical behaviour of the saddle-point equations 



(47-53). In figures || and [I| we show similar relaxation results for T = 1. As expected, 
at higher temperatures the two distributions Dj(±l, h) acquire a shape which becomes 
more like a Gaussian one, whereas in the low temperature regime the deviations from 
a Gaussian shape become important. 

To emphasise the increase in accurateness obtained by the present advanced version 
of our theory, as opposed to the simple two-parameter theory of [0, we show in figure [5] 
the simulation data and the predictions of the two versions of our theory (simple versus 
advanced) corresponding to a relaxation from a random initial state (with mo = 0), for 
T = Jq = 0. The failure of the two-parameter theory to account for the typical slowing 
down of the dynamics appears to have been amended convincingly by choosing as the 
dynamic object the full distribution Df(a,h), rather than just the magnetisation and 
the energy per spin. Since the solution of our diffusion equation (|i~7|), as depicted in 
figure ^, is obtained within the RS ansatz, this slowing down of the dynamics is not 
caused by replica symmetry breaking. 
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Figure 1: Evolution at T = of the magnetisation m and the energy per spin E, for 
Jo = (left) and Jo = 1 (right) . Solid lines: numerical simulations with N = 8000; 
dotted lines: result of solving the RS diffusion equation. 
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Figure 2: Evolution at T = of the two field distribitions Dt(a,h), for Jo = (left) 
and Jo = 1 (right) . Histograms: numerical simulations with N = 8000; lines: result of 
solving the RS diffusion equation. 
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Figure 3: Evolution at T = 1 of the magnetisation m and the energy per spin E, for 
Jo = (left) and Jo = 1 (right) . Solid lines: numerical simulations with N = 8000; 
dotted lines: result of solving the RS diffusion equation. 



J =0, rn =0.3 



J =l, rn =0.3 




Figure 4: Evolution at T = 1 of the two field distribitions Dt(a,h), for Jo = (left) 
and Jo = 1 (right) . Histograms: numerical simulations with N = 8000; lines: result of 
solving the RS diffusion equation. 
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E 




Figure 5: Comparison of simulations (N=8000, solid line), the simple two-parameter 
theory (RS stable: dotted line, RS unstable: dashed line) and the advanced theory 
(solid line), for T = J = 0. Note that the two solid lines are almost on top of each 
other at the scale shown. 



5.2 Relaxation Near the SG Transition 

One way in which we can complement the short-time results presented so far, whilst 



avoiding having to solve the saddle-point problem ( [47H53| ) for large times, is to consider 
the dynamics in the q = (paramagnetic) region. This allows us to investigate the 
relaxation near Jo = 0, T = 1 (the critical point which marks the P — ► SG transition). 
In the paramagnetic region the RS saddle-point problem can be solved, 

m = R 1 = 2J 2 R = (aH) D J A Q 2 = ((H) 

II xy 

and the diffusion equation can be expressed entirely in terms of (averages over) the 
distribution Dt(s,h) itself. Upon also making use of the invariance of the problem 
with respect to an overall spin sign change, we can write -Dt(?, h) in terms of a single 
function, the symmetric part of which is proportional to the local field distribution: 

D t is, h) = ^F t (h) (f(*H)) D = Jdy F(y)f(y) = (f(y)) F 



In terms of Ft the diffusion equation fll7|) becomes: 

J^-RW = ^[l+tanht/SiJl^t-^J-ill-tanhf/JiJlFtfiJ + ^ll-ftanhtftDr.l^Fttx) 



x[l- (tanh(f3y)) Ft ] + (y) Ft (tanh(/3y)) Ft - (y tanh(/3y)) Ft 



(66) 
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Figure 6: Relaxation of the energy per spin E for Jo = 9 = and T G {1-0, 1.1, 1.2}. 
Dots: numerical simulations with iV = 3200; lines: results of solving the RS diffusion 
equation. 



A randomly drawn initial state corresponds to 



Fq(x) = " 2 7 



JV27T 



Since ( pq ) is relatively easy to iterate numerically, we can now compare the theoretical 
predictions with the numerical data over much larger time-scales. In figure |6| we com- 
pare the result of solving ( |66| ) with numerical simulations, for t £ [0,500], in terms of 
the energy per spin E = —^(y)p t . We observe again a satisfactory agreement between 
theory and experiment. 



6 Discussion 

The present paper is the second of a sequel in which we systematically develop a 
dynamical replica theory to describe the evolution of macroscopic observables in the 
Sherrington-Kirkpatrick Q spin- glass. Our procedure for obtaining closed macroscopic 
flow equations is based on two assumptions: (i) the flow equations are self-averaging 
with respect to the realisation of the disorder, at any time, and (ii) we may assume 
equipartitioning of probability in the macroscopic sub-shells of the ensemble. The 
procedure can be shown to be exact, if the set of macroscopic observables to which it 
is applied indeed obeys closed dynamic equations. The resulting closed flow equations 
involve a saddle-point problem, to be solved at each instance of time, formulated in the 
replica language. 
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In our previous paper Q the closure procedure was applied to the observables m 
and E (the magnetisation and the energy per spin), resulting in a two-parameter dy- 
namical theory. Here we have shown how the same procedure can be succesfully applied 
to the joint spin-field distribution D(s,h), resulting in a dynamical theory describing 
an infinite number of macroscopic order parameters. The present, advanced, version 
of our theory is again by construction exact for short times, in equilibrium, and in 
the limit where the disorder is removed. Furthermore, since the joint spin-field field 
distribution specifies the underlying microscopic states in much more detail than would 
be the case by specifying only the energy and the magnetisation (i.e. more microscopic 
memory effects are being taken into acount), the equipartitioning assumption has be- 
come much weaker. We have restricted our analysis of the saddle-point equations by 
making the replica-symmetric (RS) ansatz. On the time-scales considered in our sim- 
ulation experiments, the agreement between advanced RS theory and experiment is 
quite satisfactory. For example, the slowing down missed by the two-parameter theory 
is now well accounted for, and the theory describes correctly the relaxation near the 
spin-glass transition. At this stage we need more efficient numerical procedures in order 
to extend the time-scales for which we can solve the equations of the theory. This would 
enable us to compare, for instance, with data such as the ones in |Io| , to investigate 
the possible existence of stationary states other than the one corresponding to thermal 
equilibrium, and to see whether the theory can describe the typical ageing phenomena 
observed in numerical simulations of similar mean- field spin-glass models [jllj . 

The next and final stage of our program will be to investigate for the Sherrington- 
Kirkpatrick spin-glass the effects of replica symmetry breaking on the dynamic equa- 
tions ||. Although technically non-trivial, it is a straightforward generalisation of the 
formalism developed so far. 

Finally, a relevant question which we have not yet been able to answer is whether our 
diffusion equation (17,f^) is exact (for infinitely large systems and on finite time-scales). 
There are several approaches to this problem, each of which we plan to investigate in 
the near future. The first approach is to apply our formalism to those disordered spin 
systems for which the dynamics has been solved by other means, like the non-symmetric 
SK model (in which each of the bonds is drawn independently and asymmetrically at 
random [HP; preliminary results of this study can be found in [14]), a toy model 
used in analysing the shortcomings of the previous two-parameter approach [JD^ , or the 
spherical spin-glass By definition, however, such exercises would not yet prove 

exactness in the case of the SK spin-glass. The second approach would be to try 
to derive a diffusion equation for the joint spin-field distribution, starting from the 
equations for correlation- and responsefunctions, as obtained from the path-integral 
formalism ||. The latter approach involves (rather complicated) closed equations for 
two functions C(t,t') and R(t,t'), with two real-valued arguments each (two times). 
The present formalism also involves two functions Df(l,h) and Dt(—l,h), with two 
real- valued arguments each (one time and one field). It is therefore quite imaginable 
that both formalisms constitute exact discriptions of the dynamics of the SK model. 
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